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Abstract 



Accretion disk instabilities are briefly reviewed. Some details are 
given to the short-wavelength thermal instabilities and the convec- 
tive instabilities. Time-dependent calculations of two-dimensional 
advection-dominated accretion flows are presented. 

1. Introduction 

Accretion disk is central to our understanding of many diverse astronomical 
objects and phenomena. In particular, it is believed that disk instabilities are 
fundamentally responsible for some of the luminosity variations observed in 
these systems. 

Various kinds of instabilities may arise in a rotating, shearing gas flow. For 
example, the MHD instability in a shearing flow (Balbus & Hawley 1991) is 
found to be able to survive and develop into nonlinear turbulence. In fact, an 
external magnetic fleld is not required to sustain the instability since magnetic 
held can be generated and maintained self-consistently through dynamo-effect 
within the flow (Brandenburg et al. 1995). The Balbus-Hawley instability is 
not considered here, instead, I model the small-scale turbulence by assuming 
that it introduces viscosity, which can be described by the a prescription 
(Shakura & Sunyaev 1973), 



where is the local sound speed and H is the scale-height of the disk. The 
disk is also assumed to be axisymmetric and non-gravitating, therefore the 
global non-axisymmetric Papaloizou-Pringle instability (Papaloizou & Pringle 
1984) and all the self-gravitating-related instabilities are not present here. 

An a-model accretion disk may experience the thermal and viscous insta- 
bilities (Piran 1978). They could be coherent global instabilities and are often 
applied to explain the dwarf nova outbursts and the X-ray transients (see re- 
lated discussion in this volume and references within). This type of instability 
can be analysed conveniently in terms of the Af(E) relation (where M is the 
mass accretion rate and E is the surface density of the disk). Here I consider 
the case of accretion disks around black holes. The steady state solutions or 
the thermal equilibria of accretion disks can be described in a unifled scheme 
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(see Chen et al. 1995). Depending on the magnitude of the viscosity parame- 
ter a, the coohng mechanisms, and the equation of states, there exist several 
different branches (see Fig. 1). On the Af(E) plane, it can be shown that, 
if the slope is negative, then the disk is viscously unstable (Lynden-Bell & 
Pringle 1974). The thermal instability can be examined by comparing the 
cooling and heating rates near each equilibrium curve. One may consider a 
perturbation of M (increase M but keep E unchanged), or effectively an in- 
crease of temperature, if the cooling rate exceeds the heating rate, then the 
gas cools down and the original thermal equilibrium is restored. The disk is 
therefore thermally stable. On the other hand, if the heating rate exceeds the 
cooling rate, then the gas heats up further and an instability occurs. 

a = l, M=10Mq , r=10rg a = 0.01, M=10Mg , r=10rg 
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Fig. 1. The thermal equilibria of accretion disks. Here Me = iirGM / kc is 
the Eddington limit (k = 0.34), is the heating rate, and Q~ is the cooling 
rate including both local radiative cooling and radial advective cooling. A: 
'advection' and stable, U: 'unstable', and C: 'local cooling' and stable. 

Thus, it is easily understood that for optically thick disks (large surface 
density), the thermal and viscous instabilities occur at the same time on the 
middle-branch of the S-shaped curve (or on the upper-branch of the upside 
down U-curve) when the disk is radiation pressure dominated and the ra- 
dial advection is not important. In this limit-cycle type disk evolution 
may be resulted (Taam & Fin 1984) or a transition to an optically thin disk 
state may occur (Fightman & Eardley 1974). For optically thin local cooling 
dominated disks, on the other hand, the viscous mode is stable but the ther- 
mal mode is always unstable. Therefore, a stable optically thin local cooling 
dominated accretion disk may not exist, it may either collapse to an optically 
thick disk or heat up catastrophically without bound. Finally but most impor- 
tantly, the advection dominated disks are both thermally and viscously stable 
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independent on the optically depth. 

The thermal and viscous instabilities are probably the most violent global 
instabilities an accretion disk around black hole may suffer. Under such insta- 
bilities, a disk will most probably be destroyed entirely, especially in the case 
of optically thin disks. Accretion disk models which do not have such type 
of instabilities include advection-dominated disks and gas pressure dominated 
optically thick disks. In the following, I discuss the less severe instabilities 
which may exist in advection-dominated accretion disks. 

2. Short-Wavelength Thermal Instability 

The above analysis for thermal instability, however, can only be applied to the 
long-wavelength perturbations. To examine the short-wavelength instability, 
a dispersion relation equation needs to be solved. The dispersion relation can 
be derived from the time-dependent disk equations by assuming the perturba- 
tion (of the steady state solution), 6X^ of each physical quantity, X, has the 
following functional form with respect to time (t) and space (r): 

6X/X oc e^'-'^\ (2) 

where k is the wavenumber of the modes, k = 27r/A, and ^ is, in general, 
complex. Specifically, r = Re((^) represents the growth or damping rate de- 
pending upon whether r is positive or negative, and io = Im((^) represents the 
frequency of the mode. 

I here consider only the thermal mode instability. In a recent work, Kato, 
Abramowicz, & Chen (1996) demonstrated, analytically, that short-wavelength 
thermal mode instability is still present in advection-dominated accretion disks. 
The unstable short-wavelength thermal mode has been obtained numerically 
in Chen & Taam (1993). However, detailed discussion has not been made 
there. 

By solving the dispersion relation of a transonic black hole accretion disk 
model (see Chen & Taam 1993), it is seen that the short-wavelength (A ^ 6H) 
thermal instability is a traveling mode with frequencies in the order of the 
local dynamical time-scale. Figure 2 shows the growth rate and the oscillation 
frequency of the thermal mode instability for two cases with M = 16Me and 
M = 160Me, where Me = iirGM/Kc is the Eddington hmit and k = 0.34. The 
disk is optically thick with the latter one fully advection dominated. Because 
of the nature of the propagation, the instability is expected to be weak globally 
and exist only in the inner regions near the black hole. Local properties of the 
short-wavelength thermal mode instability suggest that its non-linear behavior 
may be similar to that of the inertial-acoustic mode instability discussed in 
Chen & Taam (1995). Preliminary results of Manmoto et al. (1995, this 
volume) shows that a perturbation indeed grows due the short-wavelength 
thermal mode instability, but then it is advected into the black hole. 
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a) M=16Mp, r = 6r 
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b) M=16Mj,, r=6rg 
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c) M=160Mj,, r = 6rg 
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d) M=160Mj,, r = 6r 
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Fig. 2. The growth rate and the oscillation frequency of the thermal mode 
(solid line, the dotted line represents the viscous mode). Note the increase of 
the thermal mode growth rate for short- wavelength (A ^ 6i7), and also the 
corresponding high frequencies (uj ~ ^Ir)- 

3. Convective Instability: Time-Dependent Simulations 

In an advection-dominated disk, the specific entropy of the flow increases in- 
wards, therefore, it is likely that the convective instability is present. For a 
static rotating fluid, the convective stability condition is the H0iland crite- 
rion which states that the fluid is stable against local, axisymmetric, adiabatic 
perturbations if and only if the following two conditions are satisfied (Tassoul 
1978): 

1 df (dT\ 



dr 



^-1 Vp-VS>0, 
dp) s 
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dz dr 




> 0. 



(3a) 



(36) 
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Here (r, z) is the radial and vertical coordinates in a cylindrical coordinates, 
£ is the specific angular momentum, p is the pressure, and S is the specific 
entropy. Note that the first terms in criterion (3a) is the epicyclic frequency 
and this criterion reduces to that of Narayan & Yi (1994) for one-dimensional 
gas pressure dominated disks. 
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Fig. 3. The momentum vector field of advection-dominated fiows. The regions 
which are convectively unstable according to H0iland criterion (3ab) are indi- 
cated as grey areas. The time is in unit of Vg/c. This model has a = 0.001 and 
the initial condition is a static torus. 

We have performed two-dimensional hydrodynamical simulations for ax- 
isymmetrical advection dominated accretion fiows. 2D numerical models of 
accretion disks around black holes have been studied in different situations. 
For example, Hawley, Smarr, & Wilson (1984) simulated the formation of gas 
pressure supported tori. Eggum, Coroniti, & Katz (1987) calculated radiative 
hydrodynamical models of optically thick accretion disks. Our detailed results 
have been reported elsewhere (Igumenshchev, Chen, & Abramowicz 1996). A 
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typical flow pattern (momentum vector fleld) is shown in Figure 3. Here, we 
have used a viscosity prescription v = [2/3)acl /(Ik^ and a = 0.001. The grey 
area represents the convective unstable regions according to the H0iland cri- 
terion (3ab). We see a very efficient mixing of matter. Further analysis shows 
that this process of meridional circulations and vortices transports angular 
momentum outwards and thus the effective a becomes larger. 

4. Summary 

Advection dominated flows subject to convective instability. By time depen- 
dent simulations, we showed that this process transports angular momentum 
outwards. However, how important does it effect the dynamics of the flow is 
yet to be examined. There were some numerical simulations which showed 
that convection transports angular momentum inwards (e.g. Stones & Balbus 
1995). The difference may be caused by the different physical conditions con- 
sidered. For example, in our simulation, an a-viscosity is assumed whereas in 
Stones & Balbus (1995) no viscosity is included. 

In advection dominated flows, there is short-wavelength thermal instability. 
This type of instability is expected to be weak and localized and it may be 
termed as "flicker-type" instability (Kato et al. 1996). High resolution time- 
dependent calculation should be the future investigation. 

I wish to thank Dr. Shoji Kato for support. 
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